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The lectures describe several cosmological effects produced by neutrinos. 
. Upper and lower cosmological limits on neutrino mass are derived. The role 

that neutrinos may play in formation of large scale structure of the universe is 
described and neutrino mass limits are presented. Effects of neutrinos on cos- 
. mological background radiation and on big bang nucleosynthesis are discussed. 

I Limits on the number of neutrino flavors and mass/mixing are given. 
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^ 1 Introduction 
o 

', Of all known particles neutrinos have the weakest interactions and the smallest pos- 
p sibly nonvanishing, mass. Thanks to these properties neutrino is the second most 
^ ■ abundant particle in the universe after photons. According to observations the 



X 



number density of photons in cosmic microwave background radiation (CMBR) is 



;h I = 412/cm . In standard cosmology the number density of cosmic neutrinos can 



be expressed through 

+ np = 3n^/ll = 112/cm^ (1) 



for any neutrino flavor (z/g, t'^, and Ur), assuming that there is an equal number of 
neutrinos and antineutrinos. 

Knowing the temperature of CMBR, = 2.728 K = 2.35 ■ 10~^ eV, one can 
calculate the temperature of cosmic neutrinos: 

= = 1.95 K = 1.68 ■ 10"^ (2) 
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which is true if the neutrino mass is much smaller than their temperature, m,^ ^ T,^. 
Otherwise the parameter Tj, does not have the meaning of temperature; up to a 
constant factor it can be understood as the inverse cosmological scale factor a{t). 
Theory predicts that the spectrum of cosmic neutrinos, even massive ones, is given 
by the almost equilibrium form: 

/. = [exp(p/r-0 + l]"' (3) 

with the dimensionless chemical potential ^ = fi/T usually assumed to be negligi- 
bly small. However one should note that in the expression above p is the neutrino 



momentum, while in the equilibrium distribution there stands energy E = ym^ + p"^. 

There is a small correction to expression of the order of {m^/Td)'^ where 
is the neutrino decoupling temperature, ~ MeV - this is the temperature when 
neutrinos stopped to interact with primeval plasma. This correction appeared be- 
cause at T > Td neutrinos were in equihbrium and their distribution depended on 
E/T. Distribution of noninteracting neutrinos should be a function of pa{t). In 



the case of instantaneous decoupling it turns into f{\J{p/TY + (mi^/T^)^), while for 
non- instantaneous decoupling the dependence on mass could be different. 

More details about cosmological neutrinos can be found e.g. in a recent review 
paper 

Neutrinos are normally assumed to possess only usual weak interactions with 
iy — A)-coupling to W and Z bosons. Correspondingly, if niy = 0, only left-handed 
neutrinos, ul, i.e. those with spin anti-parallel to their momentum (and parallel for z/), 
possess this interaction, while right-handed neutrinos, un, are sterile. If rriu ^ right- 
handed neutrinos would be also coupled to intermediate bosons with the strength sup- 
pressed as {rriy/ EY and their cosmological number density would be always negligible 
since their mass is bounded from above by a few eV see eqs. ( [T3|JT^ ). If however neu- 



trinos are mixed and massive (possibly with Majorana and Dirac masses) additional 



three sterile neutrinos could be abundantly produced in the early universe 

It is known from experiment that there are at least three neutrino families (or 
flavors), Vei i^ti, and Ur- From LEP data the number of light neutrino flavors with 
rriu < mz/2 is indeed three: 

AT^ = 2.993 ±0.011 (4) 

One can find references to original experimental papers in the Review of Particle 
Physics §. 

Direct experiment limits on neutrino masses are [P] 

m^^ < 3eV, m^^ < 190 keV, m^^ < 18.2 MeV (5) 

As we will see below, cosmology allows to derive an upper limit on masses of all 
neutrino flavors similar to that presented above for m^^. 

There is a strong evidence in favor of neutrino oscillations. The best fit solutions 
to the observed neutrino anomalies indicates maximum mixing between i/^ and Ur 
with mass difference about 3 x 10^'^ eV^ (for explanation of atmospheric anomaly) 
and also large mixing between z/g and another active neutrino with mass difference 
between 10^'^ — 10^^ eV^ (for explanation of the deficit of solar neutrinos). If mass 
differences are indeed so small then masses of all active neutrinos should be below 3 eV 
and right-handed neutrinos would not be practically produced in particle interactions 
in the standard theory, but as we noted below, they may be produced by oscillations. 

Except for the above mentioned anomalies, and possibly LSND, neutrinos are 
well described by the standard electroweak theory. For a recent review of neutrino 
anomalies see e.g. ref. 

In what follows we discuss the bounds in neutrino masses that can be derived 
from the magnitude of cosmic energy density and large scale structure of the universe 
(sec. |l]). Relation between cosmological neutrinos and CMBR is considered in sec. 
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In section |^ we describe the role played by neutrinos in big bang nucleosynthesis 
(BBN) and present the limits on the number of neutrino species and possible neutrino 
degeneracy. Cosmological impact of neutrino oscillations is considered in sec. ^ The 
body of the lectures is preceded by a brief presentation of basic cosmological facts 
and essential observational data (section 0). These lectures present a shorter version 
of the recent review paper where one can find details and a long list or references, 
many of which are omitted here because of lack of space and time. 

2 A little about cosmology 

The universe is known to expand according to the Hubble law: 

V = Hr (6) 

where V is the velocity of a distant object, r is the distance to it and H = a/a is 
the Hubble constant (or better to say, Hubble parameter, since it is not constant in 
time). The present day value of H is given by 

H = 100h km/sec /Mpc (7) 

with h = 0.7 ± 0.1. There are still indications for smaller and larger values of H but 
we will not go into details here. One can find discussion of determination and values 
of this and other cosmological parameters e.g. in recent papers [Q. 

The critical or closure energy density is proportional to and is equal to: 

Pe = = 10.6 keV/cm^ (8) 

Svr 

where nipi = 1.221 GeV is the Planck mass. Contributions of different forms of 
matter into cosmological energy density is usually presented in terms of dimensionless 
parameter Qj = pj/pc- According to the data the dominant part of cosmological 
energy density is given either by vacuum energy or by an unknown form of matter 



which has negative pressure, p < — p/3, and induces an accelerated expansion (anti- 
gravity) at the present epoch. Its energy density is Qmc ~ 0.7. The total energy 
density is close to the critical value so Qtot ~ 1- The contribution of the usual 
baryonic matter, as determined from CMBR, is roughly Qbh"^ = 0.022. This is is 
consistent with determination of Qb from BBN. The remaining 0.25 is believed to 
be contributed by some unknown elementary particles (though, say, black holes are 
not excluded) weakly interacting with photons - that's why they are called dark or 
invisible matter. 

We are interested in a rather late period of the universe evolution when the tem- 
perature was in MeV range or below down to the present time. For more details 
about cosmology one can see any textbook or e.g. the recent reviews [|l], ^. Initially 
all the particles in primary plasma, photons, e"'"e~-pairs, three flavors of neutrinos 
and antineutrinos, and a little baryons were in strong thermal contact and hence had 
equilibrium distributions (P) for fermions and similar expressions with minus sign in 
front of 1 for bosons. Particle energy E should stand in this equation instead of p 
but since majority of particles are relativistic, this difference is not important. The 
energy density of massless particles in thermal equilibrium is given by the expression: 



where = 10.75 includes contribution from all mentioned above particles except for 
neglected baryons. 

At that stage the energy density was almost precisely equal to the critical one, 
with accuracy better than 10~^^, the particles were relativistic (for T > rrif.) with 
equation of state p = p/3, where p is the pressure density, and thus: 



Cross-section of neutrino interactions behaves as a^y ~ G'jpE'^ and the reaction 
rate is n/n ~ a^n ~ G%T^. Here n ~ is the particle number density and Gp = 



(9) 




(10) 
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1.166 • 10~^ GeV^ is the Fermi coupling constant. Comparing reaction rate with 
the Hubble expansion rate we can conclude that at temperatures above a few MeV 
neutrinos should be in thermal equilibrium. To be more precise one should consider 
kinetic equation governing neutrino distribution in cosmological background: 

{dt + Hpdp) Mt,p) = Hxd^fy{x, y) = hoii (11) 

where p is the neutrino momentum, x = l/a{t), y = pa/niQ, a(t) is the cosmological 
scale factor, and mo is the normalization mass which we take as itlq = 1 MeV. On 
relativistic stage when T ~ 1/a is convenient to take x — rrio/T and y = p/T. 

For an estimate of neutrino decoupling temperature we neglect inverse reactions in 
the coUision integral and assume Boltzmann statistics. Each of the neglected effects 
would enlarge the decoupling temperature by about 10-15%. In this approximation 
the kinetic equation becomes 

where D is a constant. Usually in the estimates of decoupling temperature one 
takes thermally average value of neutrino momentum, {y) = 3. If we include all 
possible reactions where neutrinos may participate then D = 80{1 + g1 + gj^) with 
ql — sin^ ^vr±l/2 and qr — sin^ 9w where sin^ 9w — 0.23 and the sign stands for 
Ve and " stands for Correspondingly the decoupling temperature determined 
with respect to the total reaction rate would be = 1.34 MeV and Tlf „ = 1.5 
MeV for Ue and u^^r respectively. 

If we take into account only (z/ — e)-interactions then D = 80{gj^ + g]^) and the 
decoupling of neutrinos from electron-positron (and photon) plasma would take place 
at T^^ = 1.87 MeV and T^^^^^ = 3.12 MeV. Above T*^ the temperatures of photon- 
electron-positron and neutrino plasma should be equal. In fact in the standard model 
equality is maintained down to T me (see below) . 
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If only z/i/-annihilation into e"'"e~ is taken into account then D = lQ{gi + 9%) and 
the decouphng temperatures would be T^^^ = 3.2 MeV and T^^^^,^ = 5.34 MeV. Below 
these temperatures the total number density of neutrinos in comoving volume could 
not change. 

When temperature dropped below electron mass, e"'"e~-pairs annihilated heating 
photons but leaving neutrinos intact. As a result of this heating become higher 
than T^, eq. and relative neutrino number density dropped with respect to pho- 
tons as given by eq. ([I|) instead of earlier existed equilibrium ratio (nj, + rip)/n^ = 3/4. 

3 Cosmological limits on neutrino mass 
3.1 Gerstein-Zeldovich limit 

Since the number density of neutrinos at the present day is known, see eq. (|l|), it is 
easy to calculate their contribution into cosmological energy density, pu = I]my^ny, 
if neutrinos are stable. Demanding that does not exceed the known value of energy 
density of matter we obtain 

m^a < 95 eV fi^/i^ ^ 14 eV (13) 

a 

where the sum is taken over all light neutrino species, a = e, p,, r. This limit was 
originally derived by Gerstein and Zeldovich in 1966. Six years later the result 
was rediscovered by Cowsik and McClelland p[. In the later paper, however, the 
photon heating by e"'"e~-annihilation was not taken into account and both helicity 
states of massive neutrinos were assumed to be equally abundant. Correspondingly 
the resulting number density of relic neutrinos was overestimated by the factor 11/2. 

If all active neutrinos are strongly mixed and their mass differences are very small 
(see the end of sec. |1]) then the limit (|13|) for an individual mass would be rrii, < 4.7 
eV. 

^More detailed discussion of such limits can be found in the recent reviews [0, ||. 
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The bound ( |13|) can be noticeably strengthened because neutrino may make only 
sub-dominant contribution to Qm- Arguments based on large scale structure forma- 
tion (see below sec. \i.'3\) lead to the conclusion that Q,y < fim/S and correspondingly: 

^m,, <5eV (14) 

a 

As noted above, in the case of small mass differences the mass bound for a single 
neutrino would be < 1.7 eV. 

3.2 Tremaine-Gunn limit 

Quantum mechanics allows to obtain a lower limit on neutrino mass if neutrinos make 
all dark matter in galaxies, especially in dwarf ones . The derivation is based on 
the fact that neutrinos are fermions and hence cannot have an arbitrary large number 
density if their energy is bounded from above to allow formation of gravitationally 
bound cluster. So to make dominant contribution into dark matter neutrino mass 
should be larger than a certain value. Gravitationally bound neutrinos would be most 
densely packed if they form degenerate gas with Fermi momentum pf = myVp- The 
Fermi velocity Vp can be determined from the virial theorem: 

= GNMgal/Rgal (15) 

where Gn = ^/^^pi is the Newton gravitational constant and Mgai and Rgai are 
respectively the mass and radius of a galaxy. 

The number density of degenerate neutrinos and equal number of antineutrinos is 
riy = p|./(37r^) and correspondingly their total mass in a galaxy is 

= 47iRl^im^n^/3 (16) 

According to observations galactic masses are dominated by invisible matter, so one 

should expect that ^ Mgai. From the equations above we find: 

, , / 300 km/sec f IkpcV^^ 
„„.80ev(^^) (^j (17) 



For dwarfs Rgai ~ 1 kpc and V ~ 100 km/sec. Correspondingly neutrinos, if they 
constitute all dark matter in such galaxies, should be rather heavy, mi, > 100 eV 
in contradiction with Gerstein-Zeldovich limit. Thus we have to conclude that dark 
matter in galaxies is dominated by some other unknown particles. 

3.3 Neutrinos and large scale structure of the universe 

Though, as we saw above, massive neutrinos cannot be dominant dark matter parti- 
cles, they may play an essential role in large scale structure formation and evolution. 
According to the accepted point of view cosmological structures have been developed 
as a result of gravitational instability of initially small primordial density perturba- 
tions. The latter presumably were generated at inflationary stage due to rising quan- 
tum fluctuations of the inflaton field. For reviews and list of references see e.g. . 



It is usually assumed that the spectrum of initial density perturbation has a simple 
power law form, i.e. Fourier transform of the density perturbations 



(5p/p),„ = J d'k5{k) (18) 
behaves as 6"^ ^ k"'. Moreover, the value of the exponent, n, is usually taken to be 1. 



It corresponds to fiat or Harrison-Zeldovich spectrum |jTT[, as indicated by inflation 
and consistent with observations. 

With the known initial perturbations and equation of state of cosmological matter 
one can calculate the shape of the evolved spectrum and to compare it with observa- 
tions. This permits to determine the properties of the cosmological dark matter. In 
the case of neutrinos density perturbations at small scales are efficiently erased as can 
be seen from the following simple arguments. Neutrinos were decoupled from plasma 
when they relativistic. The decoupling temperature is T'^ ~ MeV, while rrii, < 10 
eV. Thus after decoupling neutrinos free streamed practically with the speed of light. 
Since the flux of neutrinos from neutrino- rich regions should be larger than that from 
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neutrino-poor regions, the inhomogeneities in neutrino distribution would smoothed 
down at the scales smaller than neutrino free path, Ifs = 2tnr- Here tnr is the cosmic 
time from beginning till the moment when neutrinos became nonrelativistic. As we 
mentioned above neutrinos propagate with the speed of light, so locally their path 
is equal just to t and factor 2 came from the expansion of the universe. The mass 
contained inside Ifs is 

Mfs= 3^^^ P = ml,tf, (19) 

where we used for the cosmological energy density the critical value (H) with the 
Hubble parameter H = l/(2t), as given by eq. ([T0|). Assuming that the universe was 
dominated by relativistic matter (photons and three neutrino flavors) till neutrino 
temperature dropped down to = and taking into account that !^ 0.7T^ 

(0) we find that the mass inside the free-streaming length is 

Mfs = O.lmli/ml ^ 10^^MQ{eV/m^Y (20) 

where Mq = 2 ■ 10'^^ g is the solar mass. This result is derived for the case of one 
neutrino much heavier than the others. It would be modified in an evident way if 
neutrinos are mass degenerate. 

In such a theory the characteristic mass of the first formed objects, Mfs, is much 
larger than the mass of large galaxies, Mgai ~ lO^^M© and dark matter with such 
property is called hot dark matter (HDM). The dark matter particles for which the 
characteristic mass is smaller than the galactic mass are called cold dark matter 
(CDM) and the intermediate case is naturally called warm (WDM). In HDM model 
of structure formation large clusters of galaxies should be formed first and smaller 
structures could be created from larger ones later by their fragmentation. However 
such process demands too much time and, moreover, the observations indicate that 



smaller structures are older. Together with Tremaine-Gunn limit discussed in sec. ^ 
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it "twice" excludes neutrinos as dominant part of dark matter in the universe. How- 
ever mixed models with comparable amount of CDM and HDM are not excluded. 
Though the mystery of cosmic conspiracy - why different particles have comparable 
contribution to i7 - becomes in this case even more pronounced: 

^vac ~ (^m = ^CDM + ^HDm) ~ (21) 

From the arguments presented here one can see that the larger is the fraction 
of neutrinos in the total mass density of the universe the smaller should be power 
in cosmic structures at small scales. This permits to strengthen the upper limit 
on neutrino mass. Especially sensitive to neutrino mass are the structures at large 
red-shift z because in neutrino dominated universe small structures should form late 
and should not exist at large z. The neutrino impact on the structure formation 



was analyzed in refs. [12| with the typical limits between 1 and 5 eV. More detailed 



discussion and more references can be found in the review According to ref. ||13 
Sloan Digital Sky Survey is potentially sensitive to rriu < 0.1 eV. 

3.4 Cosmological limit on heavy neutrino mass 

If there exists fourth lepton generation then the corresponding neutrino should be 
heavier than 'mz/2 = 45 GeV to surpass the LEP result (^). If these heavy neutrinos 
are stable on cosmological time scale, >tu 10^° years, then their mass density 
may be cosmologically noticeable. Since such neutrinos are assumed to be very heavy 
their number density at decoupling should be Boltzmann suppressed and they may 
escape Gerstein-Zeldovich limit. First calculations of cosmological number density 
of massive particles were performed by Zeldovich in 1965 However his result 



contained a numerical error later corrected in ref. |T5[. The same approach was 
applied to the calculations of the number/energy density of relic heavy neutrinos 



practically simultaneously in the papers ||16| where it was found that the mass of 
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heavy neutrino should be above 2.5 GeV to be cosmologically safe. 

The number density of massive particle (neutrinos) s which survived annihilation 
is inversely proportional to the annihilation cross-section and is approximately 
given by the expression 

riy/n^ = {aannVmympiY^ (22) 

where v is the cm. velocity of the annihilating particles and is the number density 
of photons in CMBR. For relatively light neutrinos, niy <^ mz (which is not realistic 
now), the annihilation cross-section is proportional to cr ~ and the energy density 
of heavy relic neutrinos drops as According to the calculations quoted above 

= pc for nil, = 2.5 GeV. 

For higher masses > m^z, the cross-section started to drop as a ~ a^/ml 
and the cosmologically allowed window above 2.5 GeV becomes closed for niy > (3-5) 



TeV [0. However for rrii, > mw a new channel of annihilation becomes open, 

lyh + jyh^ w+w- (23) 

with the cross-section rising as m^^ The rise of the cross-section is related to 

the rise of the Yukawa couplings of Higgs boson which is necessary to ensure a large 
mass of Ufi- Correspondingly the excluded region above a few TeV becomes open 
again. However the annihilation ( p3D proceeds only in one lowest partial wave and 
the cross-section is restricted by the unitarity limit ]T^ , 

aj<7T{2J + l)/p\ (24) 

If one assumes that this limit is saturated then the large values m^^ about 100 TeV 
would be forbidden. In reality the limit should be somewhat more restrictive because 
it is natural to expect that the cross-section started to drop with rising mass of 
neutrino before it reaches the unitarity bound. However it is very difficult, if possible 
at all, to make any accurate calculations in this strong interaction regime. 
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To summarize this discussion, the cosmic energy density, p^,^, of heavy neutrinos 
with the usual weak interaction is sketched in fig. (|I]). In the region of very small 
masses the ratio of number densities Uy^jn^ does not depend upon the neutrino mass 
and linearly rises with mass. For larger masses Oann ~ "^^^ and p^^^ ~ . This 

formally opens a window for vtiy^ above 2.5 GeV. A very deep minimum in p^^ near 
m^^ = is related to the resonance enhanced cross-section around Z-pole. Above 

Z-pole the cross-section of zz/iZZ/^-annihilation into light fermions goes down with mass 
as /m^^^ (as in any normal weakly coupled gauge theory). The corresponding rise in 
p^f^ is shown by a dashed line. However for m,^^ > m\Y the contribution of the channel 
Vh^h ~^ W^W~ leads to the rise of the cross-section with increasing neutrino mass as 
c"arm ~ oP'T^t^ I f^w ■ ^his would allow keeping p^^ well below p^ for all masses above 
2.5 GeV. The behavior of pi^^, with this effect of rising cross-section included, is shown 
by the solid line up to m,^^ = 1.5 TeV. Above that value it continues as a dashed line. 
This rise with mass would break unitarity limit for partial wave amplitude when 771^^ 
reaches 1.5 TeV (or 3 TeV for Majorana neutrino). If one takes the maximum value 
of the S-wave cross-section permitted by unitarity (p4|), which scales as 1/m^^, this 
would give rise to p^^ ~ m^^ and it crosses pc at 200 TeV. This behavior is 

continued by the solid line above 1.5 TeV. However for niy^ > a few TeV the Yukawa 
coupling of to the Higgs field becomes strong and no reliable calculations of the 
annihilation cross-section has been done in this limit. Presumably the cross-section 
is much smaller than the perturbative result and the cosmological bound for mj,^ is 
close to several TeV. This possible, though not certain, behavior is presented by the 
dashed-dotted line. One should keep in mind, however, that the presented results for 
the energy density could only be true if the temperature of the universe at an early 
stage was higher than the heavy lepton mass. 
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Figure 1: Cosmological energy density of massive neutrinos Q = Pu^/Pc as a function 
of their mass measured in eV. The meaning of different hues is explained in the text. 



4 Neutrinos and CMBR 

Measurements of the angular fluctuations of CMBR, which are in continuous progress 
now, also permit to obtain valuable information about cosmic neutrinos. The spec- 
trum of fluctuations is presented in terms of Ci, the squares of the amplitudes in the 
decomposition of the temperature fluctuations in terms of spherical harmonics: 

AT 

— = Y,airnYimiOA) (25) 

Lm 



and 



= ^ E (26) 



A typical spectrum of fluctuations is presented in fig. 0(a) [20 . 



14 




Figure 2: (a) An example of an angular spectrum of CMB anisotropies with varying 
number of neutrino species, = 2,3,4. (b) The ratio of Ci for = 2, 4 relative to 
k^ = 3 (from ref. 0]) 

For low / the amplitudes Ci are practically /-independent if the spectrum of initial 
density perturbations is flat. At / ^ 200 it has a pronounced peak and a few weaker 
peaks at larger /. At / > 10^ the fluctuations are strongly damped. A detailed 
explanation of these spectral features can be found e.g. in the review These 
peaks were produced by sound waves at the earlier stage, roughly speaking at the 
moment of hydrogen recombination at T ^ 3000 K. After this moment the universe 
became transparent to CMB photons and the features existed at the moment of 
recombination became "frozen" and are observed now in the sky. The first peak of 
the largest amplitude corresponds to the "last" sound wave with the wave length 
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equal to the horizon size at recombination divided by the speed of sound which for 
relativistic plasma is = 1/ y/S: 

Ai = 4rVv^ (27) 

The other peaks correspond to maximum compression or rarefaction at the same 
moment and their phase is larger by rnr. Their amplitude is typically smaller because 
corresponding waves was generated earlier and had more time to decay. 

The physical scale l^^^^^ depends upon the expansion regime and, in particular, 
upon the fraction of relativistic matter. Because of that the peak position is sensitive 
to number of neutrino families and to neutrino mass. However this effect is rather 
weak and to a much larger extent the position of the peak is determined by the 
geometry of the universe because the angle at which we see l^^^^^ at the present 
time depends upon the curvature of space. The data strongly support spatially flat 
universe, Qtot = 1 (for an analysis of previous and new data see e.g. ref. p2|). 

More sensitive to the contribution of relativistic matter are heights of the peaks. 
The point is that at non-relativistic or, in other words, at matter-dominated (MD) 
stage gravitational potential of perturbations, ^, remains constant. Indeed the po- 
tential satisfies the Laplace equation: 

(l/a^) Ai/j ~ 6p (28) 

where a is the cosmological scale factor and 6p is the density contrast. It is known that 
density perturbations at MD-stage rise as Sp/p ~ a, while the total energy density 
decreases as p ~ 1/a'^. Correspondingly 

~ a^p {5p/ p) = const (29) 

If the cosmological expansion is not exactly non-relativistic due to presence of some 
relativistic matter (neutrinos) the gravitational potential would be time depending, 
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ip = il){t), and the sound waves (which are induced by gravity) would be amphfied - 
the effect is analogous to parametric resonance amplification. 

The heights of the peaks also depend upon the fraction of baryonic matter, VL^K^ 
since the latter makes the main contribution to the mass of acoustic oscillators, while 
dark matter particles do not experience any pressure from photons of CMB. The 
degeneracy of the impact of different cosmological parameters on the angular spec- 
trum of CMBR makes determination of these parameters much more difficult and one 
needs to invoke additional information from other pieces of astronomical data and/or 
to wait till more precise measurement of all Ci in forthcoming experiments. 

At the present day the accuracy of determination of the properties of cosmic 
neutrinos from CMBR is not very good. In ref. the upper limit on the number of 
neutrino families was found < 17 (95% confidence level) for the Hubble parameter 
h = 0.72 ± 0.08 and Qbh"^ = 0.020 ± 0.002. If these parameters are larger a larger 
fraction of relativistic energy density would be allowed and more neutrino flavors 
or other relativistic particles may exist. An additional account of the data on the 
large scale structure permitted to arrive to an interesting lower limit, > 1.5. 
Thus an independent indication (in addition to BBN) of non-vanishing cosmological 
background of massless or very light neutrinos is obtained. Of course, these results 
are not competitive with BBN (see sec. §) at the present time. However they can be 
such in near future. 

If neutrinos are massive and contribute into hot component of dark matter, their 



presence can be traced through CMBR |2J]. Both effects mentioned above, a shift of 
the peak positions and a change of their heights, manifest themselves depending on 
the fraction of hot dark matter Qhdm- Moreover the angular spectrum of CMBR is 
sensitive also to the value of neutrino mass because the latter shifts teq, the moment 
of the transition from radiation dominance to matter dominance. According to the 



paper p3 the amplitude of angular fluctuations of CMBR is 5-10% larger for 400 < 
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/ < 1000 in the mixed hot-cold dark matter (HCDM) model with Vt^ = 0.2 - 0.3 in 



comparison with the pure CDM model. An analysis of the data of ref. [^] on CMBR 
angular spectrum was performed in the paper. and the best-fit range of neutrino 
mass was found: 

= 0.04 - 2.2 eV. (30) 

An interesting effect which is related to neutrino physics in the early universe 
when the latter was about 1 sec old and the temperature was in MeV range may be 
observed in the forthcoming Planck mission if the expected accuracy at per cent level 
is achieved. It is usually assumed that at that time neutrinos had an equilibrium 
spectrum with the temperature which was initially equal to the temperature of pho- 
tons, electrons, and positrons, while somewhat later at T < rrie neutrino temperature 
dropped with respect to the photon one because e"'"e~-annihilation heated photons 
but not neutrinos since neutrinos were already decoupled from electrons and positrons 
(see eq. (||) and discussion at the end of sec. ||). However the decoupling was not 
an instantaneous process and some residual interactions between and neutrinos 
still existed at smaller temperatures. The annihilation of the hotter electron-positron 
pairs, e~^e~ —>■ uu, would heat up the neutrino component of the plasma and distort 
the neutrino spectrum. The average neutrino heating under assumption that their 
spectrum maintained equilibrium was estimated in ref. However, the approxi- 



mation of the equilibrium spectrum is significantly violated and this assumption was 
abolished in subsequent works. In the earlier papers p8| , p9| kinetic equations were 



approximately solved in Boltzmann approximation. In ref. p9[ the effect was calcu- 



lated numerically, while in ref. an approximate analytical expression was derived. 
After correction of the numerical factor 1/2 the calculated spectral distortion has the 
form: 

6L 4 E niE 



, ^ 3 • 10"^ — — - - 3 (31) 



Here Sf = f — f^^'^\ The distortion of the spectra of z/^ and Ur is approximately twice 
weaker. 

An exact numerical treatment of the problem (i.e. numerical solution of the 
integro-differential kinetic equations without any simplifying approximations) was 
conducted in the papers |^-|]32[. The accuracy of the calculations achieved in ref. 



was the highest and some difference with the results of two other papers can be pre- 
scribed to a smaller number of grids in the collision integral or to non-optimal 
distribution of them |^ (see discussion of different methods of calculations in ref. 
Recently calculations of the distortion of neutrino spectrum were done in ref. ||3^ in 
a completely different way using expansion in interpolating polynomials in momen- 
tum . The results of this work perfectly agree with those of ref. . 



One would expect that the distortion of neutrino spectrum at a per cent level would 
result in a similar distortion in the primordial abundances of light elements. However, 
this does not occur because an excess of energetic neutrinos over the equilibrium 
spectrum which would give rise to a larger neutron-to-proton ratio and to a larger 
mass fraction of primordial ^ife is compensated by an increase of the total energy 
density of Ug which acts in the opposite direction diminishing the neutron-proton 
freezing temperature and thus diminishing the n/p-iatio (see below sec. The 
net result of this distortion on '^He is at the level of 10~^. Still the observation of 
this small deviation of neutrinos from equilibrium is not impossible. The corrections 
discussed here and electromagnetic corrections of ref. could be interpreted as a 



change of N^, from 3 to 3.04. Planck mission may detect this effect but the concrete 
features depend upon the ratio of neutrino mass to the recombination temperature 



Tree = 3000 K = 0.26 eV. 



5 Neutrinos and BBN 
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5.1 BBN - brief description 

Physical processes essential for BBN took place when the temperature of the primeval 
plasma was in the interval from a few MeV down to 60-70 keV. According to eq. (|10D 
the characteristic time was respectively between 0.1 sec up to 200 sec. During this 
period light elements , ^He, "^He, and '^Li were synthesized. As one can see from 
eq. (p!OD the universe cooling rate i.e. time-temperature relation depends on the 
number of particle species in the primeval plasma, g^:. At BBN epoch the latter is 
usually parameterized as 

= 10.75 + (7/4)(iV,-3) (32) 

where the effective number of additional neutrinos AA*",^ = N^—S describes any form of 
energy present during BBN. This parameterization is precise if a non-standard energy 
has relativistic equation of state, p = p/3. In any other case AA^ becomes a function 
of time and, moreover, the impact of additional energy on primordial abundances of 
different light elements could be different from that created by the equivalent number 
of massless neutrinos. 

Building blocks for for creation of light elements were prepared in the weak inter- 
action reactions: 

n + z/e ^ p + e~, (33) 
n + ^ p + u (34) 

At high temperatures, T > 0.7 MeV, these reactions were fast in comparison with 
the universe expansion rate H and the neutron-proton ratio followed the equilibrium 
curve: 

{n/p) = exp (-Am/T) exp (-^,.J (35) 

where Am = 1.3 MeV is the neutron-proton mass difference and C,ue = f^ue/T is di- 
mensionless chemical potential of electronic neutrinos. At smaller T reactions (^^4|) 
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became effectively frozen and the ratio n/p would be constant if not the slow neutron 
decay with the life-time r„ = 885.7 ± 0.8 sec ^j. 

The temperature of the freezing of the reactions ( |5B| , |5^ ) , T„p, is determined by 



the competition of the reaction rate, F ~ G^pT^ and the expansion rate, H ~ y^T^. 
Hence T„p ~ g*^^^ ■ Larger is g^^ larger is T„p, and more neutrons would remain for 
creation of hght elements. On the other hand, the nucleosynthesis temperature, T/v5, 
does not depend upon but the time when T/yg is reached depends upon it. Larger is 
g^, shorter is this time. Correspondingly less neutrons would decay and more helium-4 
and deuterium would be created. Thus a variation of g^ acts in the same direction in 
both phenomena. In particular, an increase of by 1 leads to an increase of ^ife 
by about 5% as one can easily check using the publicly available BBN code |^ . 



Light elements begin to form from the primeval protons and neutrons when the 
temperature reached the value 

0.064 MeV 



■NS 



(36) 



1 - 0.029 In r/io 

where r/io is the ratio of baryon and photon number densities in units 10^°. A very 
small baryon-to-photon ratio makes nucleosynthesis temperature much lower than 
typical nuclear binding energies. A few years ago r] was determined from BBN itself 
and now the measurements of the angular spectrum of CMBR (see sec. ^) permit to 
find it independently. According to the latter r^io ~ 5. This result is in a reasonable 
agreement with determination of ?7io from BBN. Once T/vs is reached light elements 
are very quickly formed and practically all neutrons, that survived to this moment, 
were binded in ^He. The mass fraction of this element was about 25%, while relative 
number of deuterium nuclei (as well as ^He) with respect to hydrogen is (a few) x 10~^. 
Amount of Li is about five orders of magnitude smaller. Primordial abundances of 
light elements as functions of t^iq are presented in fig. ^. 

A slight increase of '^He as a function of 7] can be mostly explained by an increase of 
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Figure 3: Abundances of light elements (by number) '^He (by mass), and Li (by 
number) as functions of baryon-to-photon ratio r^io = lO^^ns/n^. 

the nucleosynthesis temperature (|36D with rising rj. Correspondingly nucleosynthesis 
started earlier and more neutrons survived decay. A strong decrease of the amount of 
the produced deuterium is explained by a larger probability for to meet a nucleon 
and to proceed to '^He. 

To summarize, we have seen that the primordial abundances depend upon: 

1. Number density of baryons, rj = UB/n^- 

2. Weak interaction rate; this is usually expressed in terms of the neutron life- 
time. It is interesting that a variation of weak interaction strength by a factor 
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of few would result either in complete absence of primordial '^He or in 100% 
dominance of the latter (no hydrogen) . In both cases stellar evolution would be 
quite different from what we observe. 

3. Cosmological energy density; non-standard contribution is usually parameter- 
ized as additional number of neutrino species AiVj,. 

4. Neutrino degeneracy; degeneracies of z/^ or are equivalent to a non-zero AA'^^, 
while degeneracy of has a much stronger (exponential) impact, see eq. (|35|), 
on the abundances because z/e directly enters the (n-p)-transformation reactions 



(33,34) and can shift the (?T,/p)-ratio in either direction. 



5.2 Role of neutrinos in BBN 

The role that neutrinos played in BBN is already clear from the previous section. The 
sensitivity of BBN to the number of neutrino families was first noticed by Hoyle and 
Tayler in 1964 Two years later a similar statement was made by Peebles 



More detailed calculations were performed by Shvartsman |^ in 1969 who explicitly 
stated that the data on light element abundances could be used to obtain a bound 
on the number of neutrino flavors. Another 8 years later Steigman, Schramm, and 
Gunn presented analysis of the effect with all light elements taken into account. 

The dependence of the produced deuterium and helium-4 on the number of neu- 
trino species for different values of the baryon number density rjio are presented in 
figs. |,|[ At the present day the conclusion about the allowed number of neutrino 
families varies from paper to paper between an optimistic constraint AA^;^ < 0.2 up 
to a more cautious one AA^^^ < 1. The review of these results can be found e.g. in 
refs. il|, i. 



If cosmological lepton asymmetry is non-vanishing then chemical potentials of 
neutrinos are non-zero and their energy density is higher than the energy density of 
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Figure 4: Mass fraction of ^ife as a function of the number of massless neutrino 
species. Different curves correspond to different values of the baryon-to-photon ratio 
r]iQ = lO^^ns/n^ = 2, 3, 4, 5, 6 in order of increasing hehum abundance. 

non-degenerate neutrinos. Additional effective number of neutrino species in this case 
is given by 

2 / ^ \ 2" 



AiVK) - T ? 2 (I) + (I) 



(37) 



where the sum is taken over all neutrino species a = e, fi, r. If one extra family of 
neutrinos is allowed by BBN then chemical potentials of and Vr are bounded by 
< 1-5. For electronic neutrinos the limit is much stronger, |^e| < 0.1 [ ^ . 
If combined variation of all chemical potentials is allowed such that by some con- 
spiracy an increase in production of light elements due to non-zero ^^^r is compensated 
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Figure 5: Deuterium-to-hydrogen by number as a function of the number of massless 
neutrino species. Notations are the same as in fig. (Q). 

by a positive then the hmits would be less restrictive. Using additional data from 



CMBR which allowed to fix r/io = 5 the authors of reference obtained: 



< 2.6, -0.1 < < 0.2 (3^ 



This bound disregards mixing between active neutrinos discussed below in sec. |6.3| . 
Account of the latter may make this limit considerably more restrictive, eq. (|64D. 

Distortion of spectrum of z/^ or Vr would have an impact on primordial abundances 
only through a change in total energy density of these neutrinos and is equivalent to 
a change of N^^. However distortion of the spectrum of Ve would have a stronger and 
non-monotonic influence on the abundances. If there is an increase of neutrinos in 
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high energy tail of the spectrum then the frozen n/p-ratio would increase. Indeed an 
excess of v^. in comparison with the equilibrium amount would lead to a more efficient 
destruction of neutrons in reaction (P5|), while the an excess of v would lead to more 
efficient production of neutrons in reaction (|3^) . If spectral distortion is symmetric 
for particles and antiparticles then the second reaction would dominate because the 
number density of protons is 6-7 times larger than the number density of neutrons. 
It is worth noting that spectral distortion is not necessary charge symmetric as e.g. 
could be in the case of resonant oscillations between active and sterile neutrinos. 

If spectrum of z/g has an additional power at low energy then the nj-p ratio would 
decrease because the reaction (0), where neutrons are created, is suppressed at low 
energies due to threshold of 1.8 MeV. 

If the spectrum is not distorted but the temperature of neutrinos is different from 
the temperature of photons then for a larger Ty^ the ri/p-ratio would be smaller 
because larger fraction of z/g would shift the neutron freezing temperature, T„p, to a 
smaller value. Thus an increase of the energy density of v^, has an opposite effect on 
BBN than an increase of the energy density of z/^ or z/^. 

6 Neutrino oscillations in the early universe 
6.1 Effects of medium 

There are serious reasons to believe that the mass eigenstates of neutrinos Vj {j = 
1,2,...) do not coincide with interaction eigenstates, z/^ (a = e, /i,r, ...). This mis- 
match leads to oscillation phenomenon. The interaction and mass eigenstates are 
expressed through each other by an unitary mixing matrix which in a simple case of 
two particle mixing (e.g. z/g and u^) has the form: 

z/g = z/i cos 9 + ^2 sin 9 

= —vx sin 9 \ V2 cos 9 (39) 
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If a certain neutrino flavor (interaction eigenstate) is produced in space point a tlien 
the probability to observe such state in space point h oscillates as cos[(a — h) dm? /2p\ 
where 6m? = — m\ is the mass difference squared of the mass eigenstates. 

Oscillations in matter are modified by neutrino effective potential Veff which 
up to energy factor coincides with refraction index of neutrinos in medium. The 
Schroedinger equation for neutrino vector wave function in medium can be written 
as 

idt^ = {7i^ + Veff)^ (40) 

where the matrix Tim is the free Hamiltonian which is diagonal in mass eigenstate 
basis, while effective potential, K//, is often (but not always, see below, sec. |6l^ ) 
diagonal in interaction eigenstate basis. 



Diagonal entries oiVeff for an active neutrino were calculated in ref. |^44|| : 



V}^'}^ = ±Cir7(")Gi.T3 + C^^i^ (41) 

where E is the neutrino energy, T is the plasma temperature, Gp is the Fermi coupling 
constant, a = 1/137 is the fine structure constant, and the signs "±" refer to anti- 



neutrinos and neutrinos respectively According to ref. the coefficients Cj are: 
Ci ~ 0.95, C| ~ 0.61 and C2'^ ~ 0.17 (for T < m^). These values are true in the 
limit of thermal equilibrium, otherwise these coefficients are some integrals from the 
distribution functions over momenta. The charge asymmetry of plasma is described 
by the coefficients r^^"^ which are equal to 

T]^""^ = 2r],^ + 7],^ + 7],^ +r],- r]n/2 (for u,) , (42) 
r/(^) = 2r],^ + 7],^ + r/,, - r/„/2 (for u,) , (43) 

and 7]^'^^ for u^- is obtained from eq. ( ^3]) by the interchange fi ^ r. The individual 
charge asymmetries, rjx, are defined as the ratio of the difference between particle- 
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antiparticle number densities to the number density of photons: 

Vx = {Nx - Nji) /N^ (44) 

The first term came from thermal averaging of the time component of the cur- 
rent with which neutrino interacts (Jt). This operator is odd with respect to charge 
conjugation and has different signs for neutrinos and antineutrinos. The second term 
appears because of non-locahty of neutrino interaction related to exchange of inter- 
mediate bosons. 

The numerical values of these two terms in the effective potential as well as the 
energy difference, 6E = 5m? /2E, which determines oscillation frequency in vacuum, 
are 

iV(«) = &^^GlT'E/a = 5.59 ■ 10-^° (^)' (^) MeV (45) 
= C'ir/(")G^T3 = 1.166- 10-2^^1 (^^^ (^)VeV (46) 

Description of oscillating neutrinos in the early universe in terms of wave function 
is not adequate because the effects of breaking of coherence by neutrino annihilation 
or non-forward scattering, as well as neutrino production are essential and one has to 
use density matrix formalism [Q, Equation for the evolution of the density matrix 



has the form: 

f^ = [nup]-t{n2,p} (48) 

where Tii is the effective Hamiltonian calculated in the first order in Gp (see eq. 
(^) and 7^2 is the imaginary part of the effective Hamiltonian in the second order 
in Gp. Square brackets mean commutator and curly ones mean anticommutator. 
The last term turns into the usual collision integral in non-oscillating case. In many 
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problems this term is approximated as — 7(p — Peq) where 7 is the effective strength 
of interactions and peq is proportional to unit matrix with the coefficient feq-, the 
latter being the equilibrium distribution function. Sometimes this approximation 
reasonably well describes realistic situation but in many practically interesting cases 
this is not so and more accurate form of the coherence breaking terms should be taken 
(see e.g. ref. [IH). 

6.2 Impact of active-sterile oscillations on BBN 

All three mentioned above effects of neutrino influence on BBN may be present in 
case of mixing of sterile and active neutrinos. 

1. Production of Vs by oscillations may be efficient enough to produce noticeable 
contribution into total cosmo logical energy density making Ni, > 3. 

2. If MSW resonance transition is possible a large lepton asymmetry may be 
developed in the sector of active neutrinos and if it happens to be z/g-asymmetry 
it may result in a drastic change in primordial abundances of light elements. 

3. Oscillations may distort energy spectrum of neutrinos because probability of 
transformation depends upon their energy. 

We will discuss here mostly the first phenomenon because of continuing contro- 
versy in the literature (see refs. |l|, |^). Discussion and references to effects of spectral 
distortion and asymmetry generation can be found in review [|I|. 

Let us assume that there is mixing only between two neutrinos, one active and 
one sterile, and the MSW resonance condition is not realized (this is so if is heavier 
than its active partner, 6m'^ > 0). The density matrix is 2 x 2 and has 4 elements 
which satisfy the following kinetic equations: 

Pss = S2SEI (49) 
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-S2SEI - J drlAell"^ [paaiPl)fiP2) - P(P'i)f{Pi) 
dT\Aann? [paa{,Pl)p{P2) - f {Pz) f {Pi) 

R = WI-{1/2)TR 

i = -WR-{l/2)TI+{l/2)s26E{paa~Pss) 



(50) 
(51) 
(52) 



where R and / are real and imaginary parts of the non-diagonal components of the 
density matrix, pas = R+il, SE = dm? /2E, $2 = sin26', C2 = cos 26*, W = C2SE + Veff 
(no-resonance condition means that W ^ 0), and F is the total interaction rate of 
neutrinos with all other particles. In Boltzmann approximation it is given by the 
expression ([T2|) with D = 80(1 + g1 + gj^). Fermi corrections are calculated in ref. . 



They diminish the results by 10-15%. Integration in eq. (0) is taken over the phase 
space of particles 2,3, and 4 that participate in the reaction 1 + 2^3-1-4. 

First we introduce new dimensionless variables, x and y, as defined after eq. ([TT|). 
After that the last two equations (|5T|j5^) can be solved analytically as: 



Pa 



dxi 



S26E {paa - Pss)i exp 



^ dxo 



xi H2X2 



{iW + T/2)r 



(53) 



2 Jo Hixi 

Here sub-indices 1 and 2 mean that the corresponding functions are taken at Xi or at 

X2. 

One can estimate that the maximum production rate of sterile neutrinos takes 



place at the temperature 

^prod = (10 



15) (3/y)^/3 (W/eV')^/^ MeV 



(54) 



The first number above is for mixing of Vs with z/g, while the second one is for mixing 
with or Vr- Thus if the neutrino mass difference is not too small the production 
of sterile neutrinos is efficient when T ^ H. Hence the integrals in eq. ( ^3] ) are 
exponentially dominated by upper limits and can be easily taken: 

S2 SE 



Pa 



R + il 



2W -iV 



{pa 



Ps 



(55) 
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This is essentially the stationary point approximation which is valid if F is sufficiently 
large. 

Substituting this result into eq. ( ^Uj ) we find 

If we neglect pss in comparison with paa and in comparison with W'^ then we obtain 
that the rate of production of sterile neutrinos is 

F, = jF,sin2 2^^ (57) 

where dm is the effective mixing angle in matter. Its definition is evident from com- 
parison of eqs. (pB| ) and (0). This result is twice smaller than the estimates used in 
earher papers |50| and repeated in the recent review . 

Correspondingly the limit on oscillations parameters becomes weaker: 

{SmlJeV) sin^2C- = 3.16 ■ 10-^(^?.(T;f„J/10.75)3(AiV.)2 (58) 
{6mlJeV') sin' 2^^" = 1-74 ■ 10-^((?.(T;f„,)/10.75)='(AiV.)2 (59) 

If AA^ is not very small, then a better approximation in the bounds above would be 
\n^{l~ AN) instead of (AA^)^. 

Possibly these bounds would be even weaker if we take into account that oscilla- 
tions conserve total number density of neutrinos and the latter can be changed only 
by the neutrino annihilation into e'*'e~-pairs which froze at rather high temperature 
(see discussion after eq. (|12|). 

In the resonance case a striking effect of a huge rise of lepton asymmetry of 
active neutrinos may take place |^ leading to asymmetry close to 1. Due to small 



fluctuations of baryon asymmetry this huge lepton asymmetry may strongly fluctuate 



at cosmologically large scales p2|. After neutrino decoupling these inhomogeneities 



would give rise to signiflcant fluxes of neutrinos. The latter in turn could create 
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local electric currents with non-zero vorticity which may be sources of seed magnetic 
fields [0]. Turbulent eddies generated by such fiows could also generate gravitational 
waves potentially observable in forthcoming LISA mission R|. 



6.3 Oscillations between active neutrinos 

Normally active neutrinos are in thermal equilibrium even at low temperatures and 
oscillations between them do not lead to any change in their distribution. If how- 
ever cosmological lepton asymmetry is non-zero and different for different neutrino 
flavors then the oscillations may change it and lead to equality of all asymmetries. 
Naively one would expect that in the case of large asymmetry the mixing angle in 
matter is very strongly suppressed and neutrino flavor transformations are absent. 
This is not the case however because of a large non-diagonal matrix elements of the 



effective potential. This was flrst noticed in ref. [55] and discussed in detail in series 
of papers A clear description these phenomena was presented recently in the 

paper 
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The kinetic equations used in the previous section can be easily modified to apply 
to this case. Let us consider for definiteness oscillations between and z/^. One has 
to take into account the self-interaction processes Ve^ii ^ ^e^ii and Vf.Vf, ^ v^v^j,- The 
refraction index is determined by the forward scattering amplitude and since and 

are considered to be different states of the same particle one has to include both 
processes when there is a with momentum pi in initial state and a z/g or with the 
same momentum in the final state. The processes of forward transformation Vf, ^ 
give non-diagonal contributions to refraction index. Such transformations always 
exist, even among non-oscillating particles, but only in the case of non- vanishing 
mixing the non-diagonal terms in the effective potential become observable. 
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Now effective Hamiltonian has the form: 



Htr^ = ^E[^l-^ ^j^i; j=^(/^o + ah) (60) 
where 5E = 5m? /2E and a are Pauh matrices. The elements of the Hamiltonian 



matrix (pOf ) are expressed through the integrals over momenta of the distribution 
functions of other leptons in the plasma and, in particular, of the elements of the 
density matrix of oscillating neutrinos themselves. The structure of these terms is 
essentially the same as those discussed above for mixing between active and scalar 
neutrinos, see eq. (|¥lD. The contribution of self- interaction of neutrinos and antineu- 
trinos also contains two terms. One originates from non-locality of weak interactions 
and is symmetric with respect to charge conjugation: 



h 



' sym 



j dyy' (P + P) • (61) 



The second is proportional to the charge asymmetry in the plasma and equals 

h- = ^ / dyy' (P - P) (62) 

An essential feature, specific for oscillations between active neutrinos, is the presence 
of non-diagonal terms in the Hamiltonian (or in refraction index). In the case of 
large lepton asymmetry in the sector of oscillating neutrinos, the asymmetric terms 
in the Hamiltonian strongly dominate and, as a result, the suppression of mixing 
angle in the medium, found for (z/^ — z/s)-oscillations, disappears. To be more precise 
initially the non-diagonal matrix elements are zero but they quickly rise, if the initial 
asymmetry is not too high, and soon become very large. 

A large contribution of lepton asymmetry into effective potential permits to solve 
kinetic equations for density matrix analytically and find |1|: 

= (63) 
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where 7 = T/SE and brackets in the denominator mean averaging over neutrino 
spectrum. 

The solution of this equation is straightforward. It shows that oscillations are 
not suppressed by matter effects in the presence of large lepton asymmetry. A de- 
tailed numerical investigation of oscillations between three active neutrinos in the 
early universe in presence of a large lepton asymmetry was carried out in the pa- 
per ||58|. Similar investigation both analytical and numerical was also performed in 



the papers |^ . 

An analysis of the impact of oscillating neutrinos on BBN was performed for 
the values of oscillation parameters favored by the solar and atmospheric neutrino 



anomalies [^. For the large mixing angle (LMA) solution flavor equilibrium is es- 
tablished in the early universe and all chemical potentials ^e.^t.r acquire equal values. 
The results of the calculations for this case are presented in fig. ^. Since for these 
values of the parameters, asymmetries in muonic and tauonic sectors are efficiently 
transformed into electronic asymmetry, the BBN bounds on chemical potentials are 
quite strong, 

\Ca\ < 0.07 (64) 

for any flavor a = e, /i, r. 

This result does not destroy the mechanisms of generation of seed magnetic fields 
and gravitational waves discussed at the end of sec. |6.2| . Indeed, the maximum value 
of chemical potential that could be generated by oscillations is 0.6 and may be factor 
2-3 smaller depending on the mixing parameters. If shared between three neutrino 
species these lower values would not contradict the bound (^). Moreover, for a large 
region of parameter space the essential rise of asymmetry takes place below neutron- 
proton freezing temperature and does not produce a strong effect on BBN (see e.g. 
numerical or analytical |^ calculations). 
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Figure 6: Evolution of neutrino chemical potentials for LMA case, ^13 = 0, and initial 
values = = and = —0.1. Solid and dotted curves are obtained with and 
without neutrino self-interactions respectively. 

For the LOW mixing angle solution the transformation of muon or tauon asym- 
metries into electronic one is not so efficient. The transformation started at T < 1 
MeV below interesting range for BBN. The results of calculations are presented in 
fig. 0. 

The results presented in these figures are valid for vanishing mixing angle di^ (in 
the standard parameterization of the 3 x 3-mixing matrix. An analysis of different 



non-zero values of 613 can be found in the paper |58 
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Figure 7: Evolution of the neutrino degeneracy parameters for LOW case and the 
initial values = = and = —0.1. Notations are the same as in fig. 
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